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ABSTRACT 

A summary of mid-infrared continuum luminosities arising from dust is given 
for very luminous galaxies, L/^ > 10^^ Lq, with 0.005 < z < 3.2 containing active 
galactic nuclei (AGN), including 115 obscured AGN and 60 unobscured (type 1) 
AGN. All sources have been observed with the Spitzer Infrared Spectrograph. 
Obscured AGN are defined as having optical depth r > 0.7 in the 9.7 fim silicate 
absorption feature (i.e. half of the continuum is absorbed) and having equivalent 
width of the 6.2 /xm polycyclic aromatic hydrocarbon (PAH) feature < 0.1 fim 
(to avoid sources with a significant starburst component). Unobscured AGN are 
defined as those that show silicate in emission. Luminosity uLiy{8 //m) for the 
most luminous obscured AGN is found to scale as (1+z)^'^ to z = 2.8. For un- 
obscured AGN, the scaling with redshift is similar, but luminosities uL^{8 fim) 
are approximately 3 times greater for the most luminous sources. Using both 
obscured and unobscured AGN having total infrared fluxes from IRAS, empirical 
relations are found between uL^lS //m) and L/^j. Combining these relations with 
the redshift scaling of luminosity, we conclude that the total infrared luminosi- 
ties for the most luminous obscured AGN, Lfji^AGN obscured) in L©, scale as log 
L J ^{AG'N obscured) = 12.3±0.25 + 2.6(±0.3)log(l+z), and for the most luminous 
unobscured AGN, scale as log L//j(AGNl) = 12.6(±0.15) + 2.6(±0.3)log(l+z). 
We previously determined that the most luminous starbursts scale as log L//j(SB) 
= 11.8±0.3 + 2.5(±0.3)log(H-z), indicating that the most luminous AGN are 
about 10 times more luminous than the most luminous starbursts. Results are 
consistent with obscured and unobscured AGN having the same total luminosi- 
ties with differences arising only from orientation, such that the obscured AGN 
are observed through very dusty clouds which extinct about 50% of the intrinsic 
luminosity at 8 fim. Extrapolations of observable f;y(24 ^m) to z = 6 are made us- 
ing evolution results for these luminous sources. Both obscured and unobscured 
AGN should be detected to z ~ 6 by Spitzer surveys with f,y(24 //m) > 0.3 mJy, 
even without luminosity evolution for z > 2.5. By contrast, the most luminous 
starbursts cannot be detected for z > 3, even if luminosity evolution continues 
beyond z = 2.5. 
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Introduction 



Understanding tlie initial formation and evolution of galaxies and active galactic nu- 
clei (AGN) within the early Universe is a major goal of observational cosmology. While 
average parameters of the Universe itself and the observed clustering of galaxies are we ll 



described by conclusions from microwave background anisotropics (e.g. iSpergel et al.ll2003l ). 
the transformation from initial anisotropics to individual primordial galaxies during the first 
few billion years of the Universe (redshift z > 2) is not understood. When did the first 
galaxies assemble? When did the first generation of stars occur? How did these enrich the 
interstellar medium with metals, molecules, and dust? When and how did the massive black 
holes of AGN develop? 



Observations available to date (e.g. [Richards et al.l l2006l : iMadau et al.l Il998l ) indicate 



that these formation processes maximized at 2 < z < 3 , both for AGN and star formation 
(hereinafter, we refer to luminous, star forming galaxies as "starbursts"). Understanding 
that maximum and tracing the decline of AGN and starbursts since that epoch provides 
fundamental insight on how these processes have evolved in the Universe. 

The capabihties of the Spitzer Space Telescope (Spitzer) provide powerful and indepen- 
dent tools to address these issues, for several reasons. The most luminous galaxies known 
(L/R > IO^^Lq) are the Ultraluminous Infrared Galaxies (ULIRGs, e.g. Soifer et al. 1986, 
Sanders et al. 1988), whose luminosity arises because obscuring dust reemits luminosity 
absorbed from shorter wavelengths. The total dust emission in the infrared from galaxies 
of all luminosities dominates the cosmic back ground luminosity , and this luminosity con- 



strains the tota.1 nurn ber of galaxies and AGN (jChary and Elbazll200ll : iLagache et al. 
Le Floc'h et"aDl2005h . 



2004 



The extreme extinction at optical and ultraviolet wavelengths makes it impossible in 
many cases to observe the primary luminosity sources for AGN and starbursts. Observing 
such objects in the infrared overcomes this severe extinction. Furthermore, the wide wave- 
length range of the the Infrared Spectrograph on Spitzer (IRS; Houck et al. 2004) makes it 
possible to trace the same diagnostic spectral features all the way from z = to z ~ 3. 

A variety of observing programs with Spitzer have been dedicated to understanding the 
optically-faint populations of dusty extragalactic sources which are revealed in mid-infrared 
surveys at 24 fim with the Multiband Imaging Photometer for Spitzer (MIPS; Rieke et al. 
2004) and in near-infrared surveys with the Spitzer Infrared Array Camera (IRAC; Fazio 
et al. 2004). A major objective has been to understand sources discovered with very large 
infrared to optical ratios, typically having f,y(24 /xm) > 1 mJy and R > 24 mag. [i.e. ft, (24 
/im)/f,(R) > 1000]. 
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Spectroscopic observations with the IRS have shown that these dusty sources primarily 
fall into two categories: sources with strong absorption by the 9.7 //m silicate absorption 



Houck et al. 


2005: 


Yan et al. 


2007) 



(e.g. 

are attributed to obscuration by dust, which heavily extincts the rest frame visible and 
ultraviolet fluxes while producing the mid-infrared continuum through reemission. Th ese 
heavily extincted sources have been termed "Dust Obscured Galaxies" (jPey et al.ll2008l ). 



Sources with strong PAH emission are attributed to starbursts, with the PAH arising in 
the photodissociation r e gion surrounding t he ionized gas produced by the young, hot stars 



[e.g. 



Genzel et al. 1998 



Brandl et al.ll2006l ). Sources with little PAH luminosity but strong 
silicate absorption are explained as active galactic nuc lei (AGN) surrounded by absorbing 
clouds of dust (e.g. IShi et al.ll2006l : llmanishi et al.ll2007l ). These clouds are identified primar- 
ily because of heavy absorption in the 9.7 /zm silicate feature. Dust associated with clouds 
in such AGN can also be seen in silicate emission, if the continuum observed from the hotter 



dust is not observed th rough cooler dust which produces the absorption (e.g. lHao et al.ll2005 
Schweitzer et al.ll2008l ). 



In a previous paper (jWeedman and Houckll2008l ). we summarized results for the most 
luminous dusty starbursts with < z < 2.6 discovered by Spitzer and observed with the IRS. 
In the present paper, we present an analogous summary for the most luminous dusty AGN. 
Because these AGN are often at z ~ 2, are highly luminous, and are spatially unresolved 
in infrared or ground-based optical observations, they are often called quasars. We use the 
term AGN because imaging of some examples with the Hubble Space Telescop e shows that 
they have a surrounding, spatially extended component (IBussmann et al.ll2009l ). 



In the present paper, we trace the evolution of both obscured and unobscured AGN from 
z = to z = 2.8 using a single spectral diagnostic [z^Ly(8 /im)], determine the comparisons 
of their total infrared luminosities Lm to i^L,y(8 /xm), and compare the AGN to analogous 
results for starbursts. We also predict the f;y(24 /im) for AGN and starbursts that would be 
observed to z = 6 based on average spectra of the most luminous examples of each class. 
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2. Sample Definition and Data Analysis 
2.1. Source Selection 

2.1.1. IRS Spectra of Obscured Sources 

Hundreds of IRS spectra of AGN are now available within the Spitzer archive, ranging 
from systematic observations of previously known, nearby AGN to discovery observations 
of sources with z > 2. In this paper, we assemble data from 20 different Spitzer programs 
to summarize a wide variety of AGN, including 115 heavily obscured sources classified by 
having strong absorption in the 9.7 fim silicate feature, and 60 unobscured sources (type 1 
AGN) as classified by silicate emission. 

Most high redshift AGN discovered with Spitzer and with redshifts measured by the IRS 
are obscured sources with significant silicate absorption. This is primarily a selection effect, 
because the presence of strong silicate absorption provides an unambiguous spectral feature 
for a redshift measurement. By contrast, AGN with no silicate absorption appear essentially 
as power law spectra and cannot be assigned a redshift based on the IRS spectrum. 

Our summary is intended to present an overview of a wide range of sources with IRS 
spectra which satisfy criteria of being AGN, either heavily obscured or unobscured. All AGN 
which have been observed with Spitzer are not included, especially lower luminosity, local 
objects. We attempt, however, to include all of the most luminous AGN with IRS spectra at 
all redshifts. The primary selection criteria are designed to locate "pure" AGN, both heavily 
obscured and unobscured, and to distinguish them from sources having their mid-infrared 
luminosity arising primarily from starbursts. 

Most of the new sources initially discovered by Spitzer and having IRS foUowup spectra 
are at z ~ 2, typically with f^(24 fim) > 1 mJy. The signal to noise (S/N) of such spectra 
is low, so it is important to have measurement criteria which are easily applied to such faint 
sources. In Figure 1, we illustrate examples of rest-frame spectra (with artificially high S/N) 
for the brightest sources which have been observed at z ~ 2.5, scaled to flux densities actually 
observed for the brightest examples published so far. (These are source 90 in Table 1 (below) 
for the most luminous absorbed A GN, source 37 in Table 2 (below) for the most luminous 



type 1 AGN, and source MIPS 506 (lYan et al.l 120071 ) for the most luminous starburst.) The 



shape of the average spec trum shown for type 1 AGN derives from the averages of the two 



most luminous quasars in IWeedman and Houckl (120091 ). which are also sources 37 and 49 in 



Table 2. The shape of the average spectrum of the absorbed AGN is the average of the three 
luminous ULIRGs discussed below in section 3.2 and shown in Figure 5. The shape of the 
average spectrum of the most luminous starbursts is the average of the four most luminous 
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starbursts in IWeedman and Houckl (120081 ) 



In our previous discussion of starbursts with such redshifts ( IWeedman and Houckll2008l ). 
we measured luminosity using the peak flux density of the 7.7 /xm PAH emission feature. As 
can be seen in Figure 1, this is the observed peak of the spectrum for starburst sources at z 
> 2. A similar peak at ~ 8 /im arises in obscured AGN for a very different reason, which 
is the heavy silicate absorption longw 'ard of this wavele ngth and the heavy absorption by 
ices shortward of this wavelength (e.g. ISpoon et al.ll2004l ). The flux density fi,(8 /xm) of this 
apparent continuum peak at ~ 8 /um is the best measure of luminosity of an absorbed AGN 
at redshifts z ~ 2. The precise rest-frame wavelength of this localized peak can range from 
about 7.7 fxm to 8 ^m , depending on the characteristics of the absorption. Although we 
refer to this as the 8 iJ,m peak, measurements that we report are the maximum flux density 
at ~ 8 iim , regardless of the precise wavelength. 

It is essential to distinguish AGN and starbursts among the faint sources discovered and 
observed with Spitzer. Comparing spectra for the absorbed AGN and the starburst in Figure 
1 indicates the importance of the PAH features at 6.2 ^m and 8.6 /im for judging whether 
the peak at ~ 8 fim is dominated by PAH 7.7 /im emission or by the dust continuum of 
an absorbed AGN. Various previous studies and summaries indicate that a "pure" starburst 
has rest frame equivalent width (EW) of the 6.2 /xm feature greater than 0.4 /xm, but that 
sources with E W(6.2 nm) < 0.1 /xm have their mid-infrared continuum dom inated by AGN 



luminosity (e.g. iBrandl et al.ll2006l : llmanishi et al.l 120071 : iFarrah et al.l 120071 ) 



There are many intermediate cases in which both starbursts and AGN contribute to the 



Desai et al. 


2007; 


Saiina et al. 


2007 



Sargsyan et al.ll2008l ). For individual objects, determining the relative AGN and starburst 
contributions is important. In our studies, however, both for starbursts as done previously 
and for the AGN in the present paper, we attempt to avoid uncertainties regarding composite 
sources by concentrating only on sources that appear either as pure starbursts (very strong 
PAH) or as pure AGN (very weak PAH). 

This is an empirical approach, designed to make best use of the numerous spectra of 
poor S/N. It also allows a straightforward discussion of the luminosities of AGN or of star- 
bursts, without introducing model-dependent uncertainties arising from different methods of 
deconvolving composite sources. The criterion which we adopt for including sources in the 
present paper as AGN is that sources have a measured value of EW(6.2 /im) < 0.1 fim , 
or that no detection is reported for the 6.2 /xm feature. 
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2.1.2. Evidence that Absorbed Sources are Obscured AGN 



The type 1 AGN included in this paper are all classified using available optical spectra, 
so there is no ambiguity regarding the presence of an AGN. The obscured sources suffer 
severe extinction in the optical, however, so optical spectra, even when available, would 
not necessarily reveal an obscured AGN. Our definition of obscured AGN states that in 
the absence of the PAH features from starbursts, the infrared spectrum arises from dust 
reradiation of luminosity initially produced by an AGN. For a source so hidden that no 
spectral diagnostics can escape, it is never possible to disprove the possibility that some 
of the luminosity arises from hidden starbursts. There are numerous reasons to believe, 
however, that sources characterised by strong silicate absorption and weak PAH emission 
show a spectrum which arises from dust heated by an AGN and which is then absorbed by 
intervening, cooler dust. 

The first source for which silicate absorption was observed was the prototype Seyfert 2 
galaxy NGC 1068, for which the continuum of heated dust arises from a region surrounding 
the AGN which is even smaller than the narrow line region (ITelesco et al.lll984l ). As Spitzer 
IRS spectra accumulated, it was found that type 2 (optically absorbed) AGN characteristi- 
cally show silicate absorption, wh ereas type 1 (optically transparent) AGN typically show 
silicate emission (IHao et al.l 120071 ). reinforcing the association between silicate absorption 
and obscured AGN. 

The most thoroughly observed obscured ULIRG is Markarian 231. This ULIRG un- 
ambiguously contains a type 1 A GN, immersed in heavy clouds that produce absorption in 
X-ray, ultraviolet, and visible (e.g. [Gallagher et al.ll2002l ). Half of the continuum is absorbed 



by the si licate absorption at 9.7 fim, so that the optical depth r of the absorption feature is 
T = 0.7 (lArmus et al.l 120071 ). In selecting "obscured AGN", we use a criterion that assures 
sufficiently significant absorption that the source can be recognized in any IRS spectrum that 
includes the rest-frame feature. Our criterion requires that r > 0.7, so Markarian 231 defines 
the lower limit of absorption for obscured AGN which we include in the present paper. 



Imanishi et al.l (120071 ) use measures of relative absorption in the 9.7 fim and 18 /im sil- 



icate features to show from radiative transfer arguments that the heavily absorbed ULIRGs 
require a concentrated source of dust heating, which arises naturally for an AGN but not for 
distributed starbursts. Coupled with the weak PAH features in such sources, they conclude 
that these absorbed ULIRGs are dominated by AGN luminosity. 

When the overall spectral energy distributions (SED) of the absorbed sources discovered 
by Spitzer are examined, there is a distinct relation between the SED and the presence of 
strong silicate absorption. The absorbed sources show power law SEDs from mid-infrared 
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through optical wavelengths (e.g. IWeedman et al.ll2006d : iPoUetta et al.ll2008l ). By contrast, 



selection of sources which show the stell ar absorption featu re at rest frame 1.6 /im invari- 
ably yields sources with PAH emission (IFarrah et al.ll2008l ). For large samples of Spitzer 
sources having photometry that determines the SED, the power law SEDs correlate with 
the pr esence of X-ray luminosity, leading to the conclusion that these power law so urces are 
AGN jBrand et al.lbood : IPoUetta et al.ll2007l : IPonlev et allboosl : IPiore et al.l[2008h . 

For all of the above reasons, we feel it is appropriate to interpret the sources with deep 
silicate absorption and weak PAH emission as having the mid-infrared continuum arising 
because of an AGN. As emphasized in section 2.1.1, we attempt to avoid the uncertainty 
of AGN-|-starburst composite sources by using only sources with with quantitatively weak 
PAH features. This produces a sample of sources which can be considered as the "purest" 
absorbed AGN, and it is this sample which is discussed below. 



2.2. Spitzer IRS Samples of AGN 



Absorbed sources are chosen from various Spitzer observing programs as those sources 
with IRS spectra which show r > 0.7 for the 9.7 /xm silicate absorption feature and EW(6.2 
/xm) < 0.1 for the 6.2 /im PAH emission feature. Emission sources are chosen from various 
programs as those sources with 9.7 /im silicate emission of any strength and EW(6.2 /im) 
< 0.1 for the PAH emission feature. We have not attempted to include a comprehensive list 
of low luminosity, local AGN observed with Spitzer because our objective is to locate the 
most luminous sources distributed over the largest possible redshift range. 

The various observing programs by Spitzer used to assemble this collection of IRS 
spectra of AGN are summarized below. The numbering for the list below of these observing 
programs is the same as the reference codes in Table 1, which contains the sample of obscured 
(absorbed) AGN, and Table 2, which contains the sample of unobscured (emission) AGN. 
Sources discovered with Spitzer M IPS 24 nm surveys de rive from the NOAO Deep Wide 
Field Survey (NDWFS) in Bootes (jjannuzi and Deylll999 l). the Spitzer Wide-area Infrared 
Extragalactic Survey (SWIRE, Lonsdale et al. 2003, 2004)0, and the Spitzer First Look 
Survey (FLS, Fadda et al. 2006). The sources in Tables 1 and 2 arise as follows: 

1. 16 absorbed and one silicate emission Ultraluminous Infrared Galaxi es (ULIRGs) 



discov ered by the Infrared As tronomical Satellite (IRAS) from 53 sources in iFarrah et al. 
fl2007[ ). Except for 9 sources in lArmus et al.l (120071 ) . spectra in this sample are not published, 



http : / /swire . ipac . caltech . cdu/swire / astronomers 
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so we extracted spect ra from the Spitzer archive for all 37 sources shown in Figure 21 
of iFarrah et al.l (120071 ) as having r > 0.7 and for all 4 sources shown as having silicate 
emission. Spectra were extracted as described below in section 2.3. From these spectra, the 
EW(6.2 /xm ) was measured to determine if the sources satisfy our criteria for pure AGN. 
16 absorption sources and two silicate emission sources did so, and the fy(8 /im) of these 
were measured from the extracted spectra. One of the emission sources is 3C 273, included 
below in reference 16. 



2. 12 absorbed IRAS ULIRGS from 48 sources in llmanishi et al.l (120071 ) . using those 
listed with r > 0.7 and EW(6.2 fim) < 0.1. The f,y{8 /im) are measured from the published 
spectra. 



3. 13 absorbed IRAS Faint Source Catalog (FSC) ULIRGs from 28 sources in lSargsyan et al. 



(120081 ) ■ using those hsted with r > 0.7 and EW(6.2 /im) < 0.1. The 1^(8 //m) are measured 
from the original spectral extractions. 



4. 4 absorbed local AGN from 97 sources in lShi et al.l (120061 ). These are the 4 AGN out 
of 97 published spectra which have r > 0.7; the f^{8 /im) are measured from the published 
spectra. 



5. 3 absorbed AGN and 11 silicate emission AGN from 60 sources in Weedman and Houck 



(120091 ) . selected from the flux limited sample with f,y(24 fim) > 10 mJy taken from the FLS 
and the NDWFS in Bootes. The fy(8 fim) are the published values. 

6. 9 absorbed AGN with new spectra described in section 2.3.1, below, and shown in 
Figure 2. Sources derive from sources with IRS spectra in the FLS and SWIRE Lockman 
fields having fjy(24 fim) > 5 mJy and i? > 20 mag; two additional sources satisfying these 
criteria are contained in references 5 and 12. 



7. 3 absorbed AGN from 16 sources in iBrand et al.l (l2008al ). selected as X-ray sources 
in the NDWFS Bootes field having f,y{24 /im) > 1 mJy and R > 22 mag. 



8. 2 absorbed AGN from 16 sources in iFarrah et al.l (120091 ). selected from the SWIRE 
Lockman field as sources with fiy(70 /im) > 15 mJy and i? > 23 mag. 



9. 3 absorbed AGN from 11 sources in iBrand et al.l (j2008bl ). selected from the Bootes 
field as sources with f^(70 /fm) > 30 mJy and R > 22 mag. 



10. 5 absorbed AGN from 31 sources in iHouck et al.l (120051 ) . selected from the Bootes 
field as sour c es wit h fjy(24 /im) > 1 mJy and i? > 24 mag. New extractions of sources in 
Houck et al.l (|2005l ) were done as described in section 2.3 below in order to determine if 
ambiguous cases are starburst spectra with PAH or AGN spectra with absorption. The final 
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sources selected are in Table 1, and the f;^(8 /im) are from the new spectral extractions. 



11. 5 absorbed AGN from 18 sources in IWeedman et al.l (l2006al ). selected as radio 



so urces in the FLS f ield wi th fi,(24 /im) > 1 mJy and i? > 24 mag. New extractions of sources 



m 



Weedman et al.l (j2006al ) were done as described in section 2.3 below in order to determine 



if ambiguous cases are starburst spectra with PAH or AGN spectra with absorption. The 
final sources selected are in Table 1, and the f^{8 fim) are from the new spectral extractions. 



12. 3 absorbed AGN from 20 sources in IWeedman et al.l (l2006cl ). selected from the 
SWIRE Lockman field as sources with f,y(2 4 /jm) > 1 mJy and i? > 24 mag. New extrac- 
tions of sources in IWeedman et al.l (l2006d ) were done as described in section 2.3 below in 
order to determine if ambiguous cases are starburst spectra with PAH or AGN spectra with 
absorption. The final sources selected are in Table 1, and the fj,(8 /im) are from the new 
spectral extractions. 



13. 7 absorbed AGN from lPoUetta et al.l (120081 ) . selected as the most luminous absorbed 
sources known and modeled with a dusty torus. Sources have f;y(24 /xm) > 1 mJy and R > 
24 mag. 



14. 22 absorbed AGN from 52 sources in lYan et all (120071 ) and ISaiina et al.l (120071 ) . 
selected from the FLS using the criteria fv(24 /im) > 0.9 mJy, z/f^(24 /im)/zyf;^(8 fim) > 3, 
and z/f,^(24 /im)/z/f,^(0.7 /tm) > 10. The f;^(8 /im) are measured from the published spectra 
for those sources listed has having r > 0.7 and EW(6.2 /im ) < 0.1, or no detection of 6.2 
/im . 



15. 8 absorbed AGN from the 16 sources in iMartinez-Sansigre et al.l (120081 ) not dupli- 
cated in any of the samples above, selected as "obscured quasars" in the FLS field based on 
infrared colors and radio luminosities. Absorption r and the fi,(8 /im) were measured from 
published spectra for sources without detected PAH. 



16. 5 silicate emission AGN (Table 2) from iHao et al.l (120051 ) . selected as luminous, 
optically-classified type 1 AGN. 



17. 13 silicate emission AGN (Table 2) from lSchweitzer et al.l (120081 ). selected as known, 
optically-classified type 1 AGN. 



18. The silicate emission AGN PG2112+059 (Table 2) from iMarkwick-Kemper et al. 



(I2OO7I ). 



19. 3 silicate emission AGN with new spectral extractions described in section 2.3, 
below. Sources derive from sources with IRS spectra in the FLS and SWIRE Lockman fields 
having fi,(24 /im) > 5 mJy and i? > 20 mag. which have optical redshifts and classifications 
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cited as type 1 AGN in the National Extragalactic Database (NED). 

20. 26 silicate emission AGN with new spectral extractions described in section 2.3.2, 
below. Sources are those which have optical redshifts and classifications as type 1 AGN, 
either as listed in NED or as determined from examination of spectra from the Sloan Digital 
Sky Survey (SDSS, Gunn et al. 1998), and have IRS spectra available within Spitzer Legacy 
Program 40539 (G. Helou, P.I.), a flux limited sample of sources with fj/(24 yum) > 5 mJy. 
(One source which also has i? > 20 mag. is included above in reference 19.) 

Co-ordinates, redshifts, measured values of fi.(8 ;Ltm), and the resulting luminosities 
iyLiy{8 fim) for these 175 AGN are given in Tables 1 and 2. Luminosities are determined for 
Ho = 71 kms-^Mpc-\ ^]M=0.27, and ^1^=0.73. 



2.3. New Observations and New Spectra 

2.3.1. New Obscured AGN 

The IRAS ULIRGs give large samples of absorbed AGN with z < 0.3, and the sources 
discovered by Spitzer with f;y(24 fim) > 1 mJy give large samples of absorbed AGN with 
z > 1.5. For redshifts 0.3 < z < 1.5, existing samples are much more limited. We under- 
took to enlarge such samples by seeking new sources at these intermediate redshifts. This 
was done by using intermediate flux and magnitude criteria to select sources. Because the 
absorbed sources have heavy extinction in the optical, and because we also wanted sources 
at intermediate redshift, we only searched for new sources with i? > 20 mag. Because we 
desire high luminosity sources, we only searched for new sources with f,y(24 /im) > 5 mJy. 

To obtain this new sample of absorbed AGN, we flrst ex amined all sources in the F LS 



fiFadda et al.ll2006f ) and in the SWIRE Lockman Hole survey flLonsdale et al.ll2003l . 120041 ) to 

locate sources with f;/(24 /im) > 5 mJy and R > 20 mag. Within the two surveys, 51 sources 
satisfy these criteria. Of these 51, IRS spectra are previously published (2) or available in 
the Spitzer archive (31) for 33 sources. We obtained new observations for an additional 14 
of these sources. We made new spectral extractions for the 14 new spectra and for the 31 
unpublished spectra in the archive, giving a total of 47 IRS spectra out of the 51 sources in 
the complete sample. 

Most spectra have observations in the IRS Short Low module in orders 1 and 2 (SLl 
and SL2) and with the Long Low module in orders 1 and 2 (LLl and LL2), described in 
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Houck et al.1 ^OoM These give low resolution spectral coverage fr om ~5 /xm to ~35 



Spectral extractions were done using the SMART analysis package (IHigdon et al.ll2004l ) with 
either vl5 or vl8 of the Spitzer basic calibrated data products. Extractions used all off- 
order observations as source background (i.e. all LL2 observations are coadded to produce 
the background for each nod of LLl), or alternative nod positions provide background sub- 
traction for SLl when SL2 was not available. To increase S/N, spectral extractions used an 
extraction window of average width (perpendicular to dispersion) of 4 pixels, compared to 
the standard extraction width of 8 pixels that applies for flux calibrations. These "narrow" 
extractions are corrected to accurate fluxes using extractions in the same fashion of the 
standard calibrating source Markarian 231. 

The resulting 45 spectra were examined to find those which satisfy the criteria of r > 
0.7 and EW(6.2 /im ) < 0.1. Of the 45, 9 sources satisfy these criteria; these are illustrated 
in Figure 2. Sources 73 and 111 are from archival program 30447 (G. Fazio, P.I.), source 
90 from archival program 20629 (L. Yan, P.I.), source 16 from archival program 40539 (G. 
Helou, P.I.), and sources 13,20,21,23, and 70 from our own new observations, program 50031. 



2.3.2. New Unobscured, Silicate Emission AGN 

IRS spectra of faint type 1 AGN having only broad and weak silicate emission features 
cannot determine accurate redshifts. For this reason, there are no samples of silicate emission 
AGN with redshift determinations only from IRS spectra, in contrast to the large numbers 
of absorbed AGN which are discovered with Spitzer and have redshifts from the IRS. To 
be confident of the redshift for an unobscured, silicate emission AGN, we require a redshift 
determined from an optical spectrum. 

In order to increase the sample of silicate emission AGN, we have extracted new IRS 
spectra of sources which already have optical redshifts and classifications as type 1 AGN, as 
cited in NED. These new sources arise from two criteria: 

1. The sample of 45 IRS spectra described above in section 2.3.1 from the FLS and SWIRE 
Lockman Field with fi/(24 /zm) > 5 mJy and i? > 20 mag. Three of these sources have 
optical redshifts and classifications as type 1 AGN and show silicate emission features. (All 
3 also have weak PAH, EW(6.2 /zm) < 0.1, to meet our criterion of negligible starburst 
luminosity). These sources are listed in Table 2 as reference 19. 



^The IRS was a collaborative venture between Cornell University and Ball Aerospace Corporation funded 
by NASA through the Jet Propulsion Laboratory and the Ames Research Center 
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2. All sources from Spitzer program 40539 (G. Helou, P.I.) which have optical redshifts and 
classifications as type 1 AGN, either as listed in NED or as determined from our examination 
of spectra available from the SDSS. This Spitzer Legacy Program is a fiux limited sample 
of 330 sources with f,^(24 nm) > 5 mJy having complete IRS low resolution spectra. We 
identify 25 sources which have SDSS spectra with a classification as type 1. There are an 
additional 13 sources with NED citations to classifications as type 1 AGN not derived from 
SDSS spectra. We extracted all 38 of these spectra in order to determine those sources which 
show silicate emission. Of the 25 SDSS type 1 sources, 17 show silicate emission features 
[all also have EW(6.2 fim ) < 0.1]. One of these, source 16 in Table 2, satisfies the criteria 
of set 1, above, and is listed under reference 19; the remaining 16 SDSS unobscured AGN 
are listed in Table 2 as reference 20. Of the 13 additional type 1 sources, 10 show silicate 
emission features [all also have EW(6.2 //m) < 0.1], and these 10 are also listed in Table 2 
under reference 20, with citations to the original redshift sources. 

Our only use of the spectra for these silicate emission sources with optical redshifts is 
to measure the rest frame uL^{8 /im). This measurement is listed in Table 2, so these new 
silicate emission spectra are not illustrated. 



3. Discussion 

3.1. Luminosity Evolution of Most Luminous Dusty AGN 

We first present a description of the redshift dependence on luminosity for all of the 
AGN, obscured and unobscured, using the data in Tables 1 and 2. This analysis is analo gous 



to our similar analysis of the most luminous starbursts (jWeedman and Houcld l2008l ). In 
both cases, the measure of luminosity is a consistent spectral feature (7.7 /zm PAH for 
starbursts or 8 /im continuum for AGN) which can be traced over all redshifts < z < 2.8 
using observations with the IRS. Using a single spectral feature allows the straightforward 
comparison of sources over this full redshift range. 



3.1.1. Luminosity Function of Obscured AGN 

For the obscured AGN characterized by silicate absorption, the distribution of dust 
continuum luminosities with redshift is shown in Figure 3. Because the IRS spectra allow 
redshift determinations for AGN at high redshift having strong silicate absorption, and 
because many IRAS ULIRGs at low redshift are obscured AGN, there is a good distribution 
of these obscured AGN over all redshifts < z < 2.8 (triangles and crosses in Figure 3). 
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We quantitatively determine the form of evolution for the most luminous obscured AGN 
within this sample. The form of this luminosity evolution as a function of redshift is obtained 
by using the most luminous source in each interval of 0.02 in log(l+z). There are 29 such 
intervals in the redshift range < z < 2.8, and all intervals except two have sources within 
the interval. Three very luminous ULIRGs at low redshifts in Figure 3 (the sources whose 
spectra are shown in Figure 5) are excluded from the fit, as discussed further below in section 
3.2. 

The linear (first order) least squares fit to the most luminous obscured AGN is given by 



\og[uK (8 fim)] = 44.87(±0.09) + 2.60(±0.27) log(l+z) for in ergs s-\ 

(1) 

The solid line in Figure 3 illustrates this fit to z = 2.5 for maximum luminosity within < 
z < 2.8 (uncertainties in this equation are the one sigma uncertainties arising from the least 
squares fit). 



3.1.2. Luminosity Evolution of Unobscured, Type 1 AGN 

The sample of unobscured, silicate emission, type 1 AGN discovered with IRS obser- 
vations of Spitzer sources is much smaller than the sample of obscured AGN because the 
IRS cannot obtain an unambiguous redshift for sources with broad or weak silicate emission 
features. The unobscured, silicate emission, type 1 AGN from Table 2 are plotted in Figure 
4, but we do not attempt a fit to their luminosity evolution because the redshift intervals 
are so sparsely sampled. The distribution of points indicates, however, that their form of 
evolution is similar to that for the obscured AGN in Figure 3. Trends are offset because of 
the systematically higher luminosities of the unobscured AGN, discussed further below in 
section 3.5. 

The distribution of the silicate emission AGN in Figure 3 matches the form of luminosity 
evolution for type 1 AGN in the Bootes field initially discovered as Spitzer 24 nm sources 



and subsequently found to have optical redshifts and type 1 classifications (IBrown et al. 



20061 ). Using a sarn p le of 183 sources with fi,(24 /zm) > 1 mJy and optical redshifts 1 < 



z < 5, iBrown et al.l (120061 ) determine that pure luminosity evolution fits the results with 
evolution of the form log L(z) = log L(z=0) + 1.15z - 0.34z^ + 0.03z^. This evolution curve 
is shown in Figure 3 as a dashed curve, normalized to the brightest unobscured AGN which 
has been observed in Bootes with the IRS. The similarity in shape of this curve to those 
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which we derived for the obscured AGN can be seen. The iBrown et aLl (120061 ) function 
maximizes at z = 2.5. Brown et al. emph asize that this evoluti on agrees with that found for 
type 1 quasars discovered optically, as in [Richards et al.l (120061 ) . 

The demonstrated similarity in evolution of both obscured and unobscured AGN is one 
of our most important results. It is one source of evidence, with more discussed below, that 
there are no fundamental, intrinsic differences among the most luminous AGN despite the 
very differing spectral signatures which divide obscured and unobscured AGN. 



3.1.3. Comparison to Evolution of Starhursts 
Another important result of the analysis above is that the evolution fact or of (1+z)^'^ for 



obscured AGN agrees with that determined for the most luminous starbursts (jWeedman and Houck 



20091 ). The form of evolution is determined the same way in both cases, taking the most 
luminuous source in each interval 0.02 of log (1+z) and extending to z = 2.6. In both cases, 
the most luminous sources were located within a variety of Spitzer spectroscopic programs, 
and a single spectral diagnostic was used to trace luminosities over all redshifts from to 
2.6. 



For starbursts, the parameter which was used in lWeedman and Houckl (120091 ) to measure 
evolution is the luminosity uL^, (7.7 /im) for the 7.7 ^m PAH emission feature associated 
with starbursts. This luminosity relates to the star formation rate (SFR) for starbursts, so 
luminosity evolution also represents evolution in the maximum SFR, which was found to 
scale as log(SFR) = 2.1(±0.3) + 2.5(±0.3) log(l+z), for SFR in M© yr'^ 

We note also that the form of evolution we derive for obscured and unobscured AGN 
and for starbursts is the same within the uncertainties as found for optically-luminous galax- 



ies. The comprehensive summary by iHopkins I (120041 ) of evolution parameters derived from 
optical and radio constraints shows evolution going as (1+z)^'^ and being almost entirely 
luminosity evolution, with ±0.5 uncertainty in the exponent. 

The sources which we use to track AGN evolution and starburst evolution are from 
completely different samples, each sample defined by the presence (starburst) or absence 
(AGN) of the PAH emission feature in the IRS spectra. Our observational result is that 
these two categories of luminous, dusty galaxies have evolved in the same way. The resulting 
similarity in evolution implies a close relationship between luminous AGN and luminous 
starbursts. What is the nature of this relationship? 



Attempts to understand the local ULIRGs resulted in the suggestion that the starburst 
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and A GN phenomena in lu minous, dusty galaxies occur in the same galaxy but at different 



times ( ISanders et al.lll988l ). In a simple summary of this scenario, a starburst is triggered by 
the inflow of a new supply of molecular gas caused by a galactic interaction or merger. Some 
of the gas reaches the central massive black hole and causes an increase in its luminosity, but 
the AGN is not seen until the radiation from the AGN dissipates the surrounding gas and 
dust, which may also quench the starburst. A similar scenario has been suggested for the 
extrem ely luminous, d ust obscured galaxies at high redshift, including those in the present 



paper (jPey et al.ll2008l ). In this case, contrasted to the local ULIRGs, the molecular material 



is collecting about a primordial galaxy in the process of formation. 

That the luminous AGN and starburst populations track each other so closely in lu- 
minosity evolution does not require, however, that either phenomena causes the other. If 
both luminous AGN and luminous starbursts arise because of galactic interactions or other 
events which supply fresh gas and dust into the gravitational potential of a galaxy, both 
can be triggered by the same event. This triggering does not require that the AGN affects 
the starburst, or that the starburst affects the AGN; they can be coincidental in the same 
galaxy because a galactic interaction stimulates both. Whatever the process, these Spitzer 
IRS results show that energizing sources for dusty ULIRGS, whether starburst or AGN, have 
declined together through the past 11 billion years, as the typical luminosity has faded by a 
factor of ~ 20. 



3.2. The Extreme Local ULIRGs 

One exception should be noted to the similarities described above among the evolution 
of both categories of AGN and of starbursts. For starbursts and unobscured (type 1) AGN, 
there are no sources with z < 0.5 having luminosities comparable to the most luminous 
sources with z > 2. There is a distinctive difference for the most luminous, obscured AGN. 
The three most luminous, obscured ULIRGs in Figure 3, all with z < 0.6, are systematically 
more luminous by a factor of ~ 3 than other local ULIRGs and come within a factor of 3 of 
the most luminous absorbed AGN found anywhere in the Universe. 

The rest-frame spectra of these three sources (all discovered by IRAS) are shown in 
Figure 5, labeled by the running numbers in Table 2. The spectra illustrate that all of these 
are deeply absorbed, so that the high luminosities z/Ly(8 fim) arise despite extinction that 
may decrease the apparent 8 /im luminosities. The weakness of the 6.2 fim PAH feature is 
evidence that the high luminosities at 8 /im do not arise from a starburst component that 
enhances the apparent 8 /im luminosity with 7.7 /im PAH emission. The spectra in Figure 5 
also show how silicate absorption > 8 fim and absorption by ices < 8 /tm produces a localized 
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continuum peak at ~ 8 fim. The average shape of these spectra shown in Figure 1 is very 
similar to the spect r a of the most lumi nous absorbed AGN at z 
PoUetta et allboosl : ISaiina et allboOTh . 



2 (e.g. iHouck et al.ll2005 



For luminous starbursts and unobscured AGN, luminosity evolution can be explained 
by a systematic change in the amount of gas that falls deep within the potential well of the 
galaxy; unobscured AGN luminosity is controlled by the fraction of this gas that approaches 
the central massive black hole, and starburst luminosity is controlled by the total amount 
of gas in the circumnuclear molecular clouds. As a consequence, the luminosity evolution of 
type 1 AGN and of starbursts both is the result of having decreasing amounts of accreting 
material at lower redshifts. 

For the luminous, obscured AGN, the presence of three low redshift sources which depart 
from the systematic trend indicates that there may be an additional factor other than the 
amount of accreted material which sometimes controls the luminosity i^L,y(8 fim) of these 
sources. It is notable, however, that the unusually high luminosities of these 3 sources when 
measured by uL^{8 /xm) are not so unusual when measured by total infrared luminosity L/^. 
The 2 IRAS ULIRGs in Figure 5 have measures of L/^ (sources 3 and 32, excluding the 
FSC source), and the ratio Lif>/uL„{8 /im) is shown in Figure 6. The 2 unusually luminous 
ULIRGs have distinctly small values for this ratio, indicating that they are more similar to 
other ULIRGs when measured in Lm. 

One possible explanation for the unusually large z/Lj,(8 /im) and unusually small ratio 
L/ij/ lyLj^^S /im) may be the location and distribution of the absorbing dust clouds. Absorbing 
clouds concentrated unusually close to the AGN could give rise to warmer dust and preferen- 
tial emission at shorter wavelengths, compared to sources with cooler, more distant clouds. 
We can consistently interpret our general results for both obscured and unobscured AGN in 
context of the simple "unified model" (section 3.5, below) containing a dusty, absorbing torus 
with comparable opening angles in all AGN. These three sources indicate, ho wever, that this 
model becomes more complex for the most hea vily absorbed sources (e.g. iLevenson et al. 



20071 : llmanishi et al.l 120071 : iPoUetta et al.ll2008l ). For this reason, we determine the form of 
luminosity evolution in section 3.1.1. for obscured sources without including these three 
anomalous sources. 



3.3. Total Infrared Luminosities of AGN and Starbursts 



We have utilized uL^lS /im) as a measure of continuum luminosity from emitting dust 
because of observational necessity; this parameter is easily measured over a wide range of 
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redshifts. Of course, a more fundamental measurement is the total dust-emitted luminosity, 
L//J, which is dominated by luminosity arising at far-infrared wavelengths from dust that 
is cooler than the dust which produces the mid-infrared luminosity at 8 fim. For those 
bright sources which have total infrared observations with IRAS, we can determine the ratio 
L/r/z/L^(8 ^m). 

For the obscured ULIRGs in Tabl e 1, the resul t is sh own in Figure 6, us i ng the uL^^S 
Hm) from Table 1 and the L/^ from iFarrah et al.l (120071 ) or llmanishi et al.l (120071 ). The 
median ratio log [Ljji/iyL^{8 /im)] = 0.95±0.25, which gives that the median correction for 
obscured ULIRGs is log Ljr = log z/Lj,(8 /im) + 0.95±0.25. (Because of the small number of 
points, we do not attempt to fit a relation having luminosity dependence.). Assuming this 
also applies to other obscured AGN in our sample, this yields that for obscured AGN, 



log L J ji^AGN obscured) = log ^^'u{^ yum) " 32.63±0.25, for Lj^ in Lq and z/Lj,(8 /im) in ergs 



-1 



(2) 

This is simply an empirical result based on local IRAS ULIRGs having observed data 
that include far-infrared wavelengths. Efforts have begun to obtain total luminosities for the 
high redshift, absorbed sources di scovered by Spitze r (including many sources in Table 1) 



using submillimeter observations (ISajina et al.l 120071 ): these authors give results in the rest 



frame compared to the broad band measurement that would be made with the 8 /xm filter 



of the Spitzer Infrared Array Camera (IRAC). The empirical result in ISajina et al.l (120071 ) 



is log LiR = (2.79± 0.36) + (0.83± 0.03)log z/L^(IRAC 8 /zm), for luminosities in L©. 

Applying the synthetic photometry tool in SMART to the average spectrum of the ab- 
sorbed source in Figure 1 indicates that the flux density observed with the IRAC 8 //m filter 
is 0.68 of the monochromatic f^{8 /im) for an absorbed source with a spectrum as in Figure 
1. Using this transformation, the result from Sajina et al. becomes log L/r = 2.65 -|- 0.83 log 
i^Lj,(8 /im). The median log i^Lj,(8 //m) in Figure 6 is 45.0 (ergs s~^) or 11.41 (Lq), which 
yields log Lift = 12.12 using the Sajina et al. transformation. Using our transformation 
in Equation 3, a source with log z/Ly(8 /im) = 45.0 (ergs s^^) would have log Lm = 12.37 
(Lq). The difference of 0.25 between these two independent determinations of log Ljr is well 
within the stated uncertainties of the two relations. This comparison gives confidence that 
the transformation between L/^ and i/L^lS fim) is similar among all of the absorbed AGN 
in Table 1 and Figure 3. 

For the unobscured, silicate emission, type 1 AGN in Table 2, we determine Ljr for 
those which have total infrared luminosities measured at all IRAS wavelengths, taken from 
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the st udy of PG quasar s in ISchweitzer et al.l (120081 ) and including 3C 273 from iHao et al. 
( 120051 ) and iFarrah et al.l (120071 ) . Source luminosities Ljpt are in Table 2 and the comparison 
to uL^lS /im) are in Figure 7. Asterisks show values for sources without measured PAH 
features and diamonds for sources with weak PAH features, although there is no systematic 
difference in the correction. The median ratio log [Lir/i>L^{8 /im)] = 0.74±0.15. This 
indicates that, for unobscured AGN, 



log L/r(AGN1) = log z/Ly(8 /im) - 32.84±0.15 for Ljr in and z/L,^(8 /im) in ergs s 



-1 



(3) 



For starbursts calibrated in iHouck et al.l (120071 ). the analogous result is log Lm 
vhyiy .7 /im) + 0.78±0.2 in ergs s~^, which gives that 



log 



log L//{(SB)= log i^Lj,(7.7 /im) - 32.80±0.2 for in and z/Ly(7.7 /i.m) in ergs s ^ . 

(4) 

The corrections to L//? appear similar for AGN and starbursts, but this is not the case if only 
the dust emission is measured. For starbursts, the /im) measures PAH emission. The 

actual dust continuum is much weaker. For a pure starburst, such as the lower spectrum in 
Figure 1, the dust continuum can be estimated by fitting the continuum on either side of the 
broad PAH c omplex, and this c ontinuum is about a factor of ten weaker than the peak at 



7.7 /im (e.g. iBrandl et al.ll2006l ). A more appropriate rest wavelength at which to compare 
intrinsic dust continua is 15 /tm, beyond the silicate absorption and not contaminated by 
PAH. But we cannot observe this rest wavelength with the IRS for sources with z > 1.2. 

Using the transformations in equations (2), (3), and (4) together with the form of evo- 
lution in equation (1) from section 3.1.1., we can determine the total luminosity L/^j as a 
function of redshift for the most luminous AGN. This can be compar ed with the result pre 



yiously derived in the same fashion for the most luminous starbursts (IWeedman and Houck 



20081 ). Results are: 



a) For the most luminous obscured AGN, 



log UR{kG^obscured) = 12.3±0.25 + 2.6(±0.3)log(l+z) for Ur in L, 



(5) 
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b) For the most luminous unobscured, type 1 AGN, 



log L/k(AGN1) = 12.6±0.15 + 2.6(±0.3)log(l+z) in L©. 

(6) 

c) For the most luminous starbursts, 

log L7r(SB) = 11.8±0.3 + 2.5(±0.3)log(l+z) in L©. 

(7) 

The results in equations (5) and (6) summarize once again the fundamental similarities 
between obscured and unobscured AGN, with their bolometric luminosities being the same 
within the uncertainties. The results also indicate that the most luminous sources in the 
Universe arise from AGN and not from starbursts. 



3.4. Luminosity Functions for the Most Luminous Dusty Galaxies at z = 2.5 

3.4- 1- Luminosity Function of Obscured AGN 



Since the discoveries of IRAS, it has been known that the most luminous galaxies have 
their luminosity domi nated by continuum radiation in the infrared, and this led to the 
definition of ULIRGs (ISoifer et al.lll986l : ISanders et al.lll988l ). The sources summarized in 
this paper represent the most luminous AGN known at any redshift, but our data do not 
extend sufficiently deep to determine luminosity functions extending through a large range 
of luminosity. 

The present sample can only determine the brightest end of the luminosity function 
at high redshift for the obscured AGN. The obscured AGN which have been found in the 
Bootes, FLS, and Lockman Hole survey areas represent a reasonably complete sample at 
z ~ 2 for sources with fiy(24 /im) > 1 mJy within a total sky area of 22 deg^. This is 
because of the many programs described in section 2.2 to discover these sources. The most 
significant selection effect depends on redshift. The peak of the spectrum at 8 //m (Figure 
1) means that the f,y(24 /im) is dependent on redshift as this peak moves through the MIPS 
24 /um filter band. This explains why so many sources have been found with z ~ 2. 



-20- 



In Figure 8, the distribution of redshifts and source luminosities for the obscured AGN 
in Table 1 is shown for 1.5 < z < 3.0. This Figure demonstrates the preferential selection 
of sources with 1.7 < z < 2.7. We use this distribution of sources to determine a luminosity 
function applying to the luminosities which are shown. 

We divide the results in Figure 8 into two redshift windows, 1.7 < z < 2.2 and 2.3 < z 
< 2.7, to make independent estimates of the luminosity function in each window. Sources 
over the luminosity range 12.9 < log L/^ < 13.9 (Lq) are counted in each interval of 0.4 
in log L/zj for each redshift window. Luminosities for sources in the lower redshift window 
centered at z = 1.95 are scaled by an evolution factor (1+z)^'^ to match luminosities in the 
higher redshift window centered at z = 2.5, so that the resulting luminosity function applies 
for z = 2.5. Total volumes contained within the two redshift intervals are determined by 
scaling to the 22 deg^ included in the surveys which discovered these sources and using the 
cosmology described in the notes to Table 1. 

The resulting luminosity function is shown in Figure 9. For comparison, we also show 
the bright end of the local luminosity function for all IRAS galaxies from Soifer, Neugebauer 
and Houck (1987), scaled to z = 2.5 using the evolution factor (1+z)^-^ which we have 
determined for the dusty AGN. The IRAS luminosity function includes all galaxies and 
makes no distinction among obscured AGN, unobscured AGN, or starbursts. From the 
Spitzer results for ULIRG spectra, in section 3.2, we can conclude that the IRAS sources 
which define the bright end arc obscured AGN. The agreement between the extrapolated 
IRAS luminosity function and our new determination from Spitzer sources is independent 
confirmation of the validity of the evolution factor, and of the validity of the empirical 
transformation from vL^^^S //m) to Ljr. 

3.4-2. Comparison to Unobscured AGN and Starbursts 

To understand the total infrared luminosity density contributed by obscured AGN, un- 
obscured AGN, and starbursts, it is necessary to understand the relative luminosity functions 
and space densities for these three categories of sources. Results on this question will be fur- 
ther refined as increased understanding develops on the classification and redshift of the large 
numbers of sources in Spitzer surveys, and how these sources relate to AGN and starbursts 
found with other techniques. 

For now, we can qualitatively scale the luminosity function in Figure 9 to estimate the 
analogous luminosity functions for unobscured AGN and starbursts just by comparing the 
fluxes and numbers of such sources found in the same surveys that reveal the obscured AGN. 



- 21 - 



Among the programs for IRS follow up of Spitzer sources with f;y(24 /ii n) > 1 mJy chosen 



without considering the near infrared spectral energy distribution s (e.g. iHouck et al.l 12005 
Yan et al.ll2007l : IWeedman et al.ll2006al : IWeedman and Houckll2009l ) , the numbers of obscured 
AGN (found with silicate absorption) and the number of unobscured AGN (featureless or 
silicate emission sources) are comparable. The total number of these AGN near z ~ 2 is 
approximately the same as the total number of starbursts classified by PAH emission. PAH 
sources dominate, however, in samples deliberately sele cted using near infrared IRAC colors 



to see k the redshifted 1.6 yum stellar absorption feature (IWeedman et al.ll2006d : iFarrah et al. 
2008h . 



More quantitative results comparing the numbers of AGN and starbursts arise by classi- 
fying the spectral shapes of sources having MIPS and IRAC fluxes. Sources with power law 



spectral energy distribu t ions consistently h ave other characteristics of AGN (jPoUetta et al. 



20081 : Brand et al.ll2006l : IPonlev et al.ll2008[ ). The AGN fraction determined i n this way de 



creas es from ~ 40% at fiy(24 //m) > 1 mJy to ~ 10% at fi,(24 yum) ~ 0.1 mJy (IDonley et al. 



20081 ) . Such surveys also in dicate that the nurn bers of obscured and unobscured AGN are 
similar, to the survey limits (iHickox et al.l 120071 ). 



We can conclude, therefore, that the space densities at z ~ 2 for the most luminous 
obscured AGN, unobscured AGN, and starbursts are similar. Their relative luminosity 
densities then scale as the relative luminosities. These relative can be estimated by 
scaling the brightest examples at z = 2.5, as in Figure 1, and using the transformations to 
Liji in equations 2-4. Determining from Figure 1 that f;^(8 /im,AGNl) : fy(8 fim,AG'N obscured) 
: f,(7.7 ^m,SB) = 6 : 2.5 : 1, the result is L,fi(AGNl) : Ljr{AGN obscured) : Ur{SB) = 5.5 
: 3.7 : 1.0. 

If relative luminosity densities scale the same as these relative luminosities, this result 
means that AGN of both categories together produce about ten times the luminosity density 
as do starbursts. We emphasize, however, that this conclusion applies only at the very highest 
luminosities and does not represent a measurement over the full luminosity functions of AGN 
and starbursts. 



3.5. Unified Model and Extinction of AGN 

The transformations from i^Lj,(8 /xm) to L/^ used above are empirical and do not apply 
any corrections for extinction of z/L,^(8 ^m). Nevertheless, understanding extinction at 8 
fim is important for understanding the geometry of the sources and the amount of dust they 
contain. Extinction can be estimated under the simple assumption of an absorbing screen. 
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using mea surements of absorpt ion r at 9.7 /im and scaling to absorption by the same dust 
at 8 //m (jPraine and Lill200ll ). There are two primary uncertainties with this approach. 
One is that absorption is not taking place in a screen, but is taking place because of complex 
radiative transfer effects within dust clouds which are hotter on one side and cooler on the 



other (e.g. iLevenson et al.ll2007l : llmanishi et al.ll2007l ). The other uncertainty arises because 
the 9.7 /im absorption is often so deep that an accurate measure of r is not possible. 

A number of studies have verified the v alidity of the " u nified model" of AGN for the 



infra r ed dust emission and absorption (e.g. IShi et al.l l2006l : iHao et al.l l2005l : iHickox et al. 



20071 : iMaiolino et al.l 120071 ). Originally conceived to explain the differences between opti- 
cal classification of type 1 AGN and type 2 AGN as arising only from orientation effects 
(lAntonucci and Milled Il985l ). the model can be applied to the infrared dust continuum. In 
this case, the observed i^Lj,(8 /im) depends on viewing angle. If absorbing clouds are viewed 
from the cool side, then L/i(;/z/Ly(8 jj,m) is measured by the value for obscured sources. If 
viewed so that the hot side of the clouds is directly observed, L/ij/z/Ly(8 /im) is measured 
by the value for unobscured sources. 

If both unobscured and obscured AGN have intrinsically the same luminosities L//j, 
then the transformations from z/Li^(8 //m) to L/^j derived above in section 3.3 differ only 
because the fhyi^ /im) suffers more extinction for the obscured sources. The total L/^j is 
independent of viewing angle because it is dominated by far infrared luminosity which is not 
subject to significant extinction regardless of viewing angle. 

Comparing the transformations from fim) to Lj^ in equations (2) and (3), the 

results show that the correction is different by a factor of 1.6 between obscured and unob- 
scured AGN, in the sense that z/L;^(8 /tm) is more luminous relative to Lm in unobscured, 
type 1 AGN compared to the obscured AGN. This result is as expected in the unified model 
if the obscured sources are viewed through dust clouds and indicates a typical extinction at 
8 /tm by a factor of 1.6 in the obscured AGN, or that ~ 60% of the intrinsic 8 /tm luminosity 
emerges from obscured sources. 

An independent estimate of the extinction at 8 /im arises by comparing the observed 
luminosities of obscured and unobscured AGN, as in Figures 1 and 4. If the intrinsic lumi- 
nosities for obscured, silicate absorbed sources and the unobscured, silicate emission sources 
are the same, the scaling difference in i/L,y{8 /im) is a measure of how much luminosity has 
been extincted at 8 /im for the absorbed sources. This scaling difference in Figures 1 and 4 
is about 2.5, i.e. the unabsorbed sources are typically two to three times as luminous at 8 
/im as the absorbed sources. This result implies that ~ 40% of the intrinsic 8 /im luminosity 
emerges from the obscured sources. 
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3.6. Tracing Dusty AGN and Starbursts to z = 6 

The presence of large amounts of dust within AGN and starbursts in the early Universe 
raises many questions. Does the luminosity evolution observed to z = 2.6 continue to higher 
redshifts? When and how was the dust formed? When was the Universe most obscured? 
Were massive black holes present before or after dust formation? Which came first, starbursts 
or AGN? Progress toward answering such questions depends on tracking the dust content of 
the Universe as a function of redshift. 

The determination of luminosity evolution discussed in previous sections extends to 
z = 2.6 only because of the limit in the spectral response of the IRS. At z = 2.6, the 8 
/im peak of the absorbed continuum is at ~ 29 /xm , and the deepest absorption is at ~ 35 
Hm . For faint sources with poor S/N, these are approximately the reddest wavelengths for 
which absorbed sources can be confidently identified spectroscopically. A crucial question, 
therefore, is whether the luminosity evolution observed to z = 2.6 continues beyond that 
redshift. Even though infrared spectra for redshift measurements cannot be obtained, we 
can ask whether such sources would exist within the Spitzer MIPS surveys at 24 /xm . This 
is done by using the results described above regarding the luminosities and evolution of these 
sources to predict the fluxes which would be observed at 24 /im as redshift increases. 

The results are shown in Figure 10. Asterisks are the most luminous unobscured, type 
1 AGN, crosses are the most luminous absorbed AGN, and diamonds are the most luminous 
starbursts, all scaled to the brightest known examples shown in Figure 1. For each set, the 
upper curve shows expected fy(24 //m) if source luminosities continue to scale as (l+z)^-^ to 
z = 6, and the lower curve shows expected ii,{24 p, m) if there is no luni inosity evolution for z 



> 2.5. (The latter is the result for type 1 AGN by lBrown et al.l (120061 ) . who determine that 
the space density of luminous type 1 sources maximizes at z = 2.5.) The horizontal line in 
Figure 9 is the MIPS limit of f,y(24 //m) =0.3 mJy for typical wide area surveys such as in 
Bootes, the FLS, and SWIRE. (Smaller area surveys have reached f;y(24 yum) < 0.1 mJy.) 

Results in Figure 10 show that even if luminosity evolution continues for z > 2.5, we 
could not expect these MIPS surveys to identify starbursts with z > 3. The rapid fading of 
starbursts arises in part because the strong PAH 7.7 fim emission feature moves longward 
of the MIPS bandpass. Spitzer studies are revealing large numbers of dusty starbursts with 
z < 3 by various survey techniques and IRS follow up, and these studies should eventually 
lead to a well quantified measure of star formation densities for z < 3. For now, tracking 
dust to earlier epochs must rely on AGN. 

Figure 10 shows that both obscured AGN and unobscured (type 1) AGN should be 
found with f,y(24 fim) > 0.3 mJy to z ~ 6, even without continued luminosity evolution for 
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z > 2.5. In fact, a number of type 1 AGN have already been discovered using f,y(24 /im) in 
the Bootes field at these redshifts (IBrown et al.ll2006l : ICool et al.ll2006l ). These are shown 
as squares in the Figure. (Redshifts are not published for the quasars in Brown et al., but 
redshifts and fluxes for sources are shown in Figure 8 of that paper, and all sources with z 
> 3 are shown in Figure 9.) That the f,y(24 fim) for all of these unobscured AGN in Figure 
10 fall below the extrapolations from the brightest source at z = 2.5 is further evidence that 
the luminosity evolution of unobscured AGN turns down for z > 2.5, as already concluded 
from optical data. 

It is important to emphasize that the evidence for this turnover of luminosity evolution 
in Figure 10 derives from the luminosity of the dust-emitting continuum. One possible 
suggestion for the maximal optimal luminosities of type 1 AGN at z ~ 2.5 when observed 
in the rest-frame ultraviolet is that dust and extinction were less at that epoch. But the 
evidence of the present analysis is that dust luminosity for unobscured AGN also shows a 
similar maximum, so the explanation for this maximum is not the absence of dust. 

It is crucial to find more high redshift examples of either category of AGN, because 
these trace dust to the most d istant epochs at w hich it can be detected. The presence of the 
two quasars in Figure 9 from ICool et al.l (120061 ) with z = 5.5 is evidence that dust existed 
when the Universe was less than 1 billion years old. 

It will be much more difficult to find obscured AGN for z > 3 than to find the unobscured, 
type 1 AGN. The obscured AGN are much fainter optically than the unobscured AGN; based 
on the R magnitudes of the obscured AGN at z ~ 2.5, we would expect R > 25. Obscured 
AGN with z > 3 may explain many of the sou rces having f^(24 um) > 1 mJy and i? > 25 mag 
which show no features in the IRS spectra (IWeedman et al.l l2006bl ). The best candidates 
for such sources are those objects with power law spectra in the MIPS and IRAC bands but 
optical magnitudes fainter than 25. The only way to determine redshifts for these will be 
with optical and near- infrared spect roscopy, but this is challenging at such magnitudes (e.g. 
Brand et ahlboO?! : lOesai et al.lboosh . 



3.7. Galaxy Formation at High Redshift 



That galaxy luminosities peak at z ~ 2.5 indicat es that this epo ch corresponds to 
a maximum for the assembly of massive galaxies (e.g. iDey et al.l l2008l ). We can use the 
numerical values for luminosities and SFRs derived above to compare with other, independent 
conclusions concerning the nature of the most luminous galaxies at such redshifts. 



The evolution of maximum SFR for starbursts to z = 2.5 reviewed in section 3.3 and the 
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absence of starbursts with greater SFR at z > 2.5 discussed in secti on 3.6 indicate that the 



maximum SFR during the epoch 2.5 < z < 6 is SFR < 3500 Mq yr'^ iBussmann et all tOQ% 
derive stellar masses for dust obscured galaxies at z ~ 2 using HST observations that resolve 
the extended stellar component of the galaxy. Bussman et al. determine a median lower limit 
to the stellar mass of 3x10^° M© (depending on uncertainties for extinction corrections) for 30 
obscured galaxies similar to the sources in Table 1. At the maximum SFR which we derive, 
such a mass could accumulate within ~ 10^ yr. This conclusion indicates that starbursts 
at the luminosities we observe need have only short duration to explain the observed stellar 
masses of galaxies at z ~ 2.5. 

We can also ask whether the necessary mass for the central massive black hole (MBH) 
presumed to power the accreting AGN is consistent with the SFRs, as required if the remains 
of stars collect within the galactic nucleus to form the accreting MBH. Many HST observa- 
tions of local galaxies indicate a consistent and systematic result that every galaxy contains 



a cen tral MBH with mass ~ 0.2% of the total virial mass of the galaxy (e.g. iFerrarese et al. 



20061 ) (these authors suggest the term "Central Massive Object", or CMO, to allow a con- 



centration of smaller black holes rather than a single MBH). 

Equations (5) and (6) indicate that the most luminous AGN at z = 2.5 has Lj^ = 10^^'^ 
Lq. If the Eddington limit applies, 10^^ ergs s~^ is produced per solar mass of the accreting 
massive black hole (MBH). In this case, the minimum accreting mass for the most luminous 
AGN at z = 2.5 is then 3x10^ Mq. This number for the high redshift, du sty AGN is the 



same as the most massive local CMO observed by HST, which is shown in IFerrarese et al. 
J2006h to have a CMO mass of 3x10^ Mq. 



Using the 0.2% scaling of Ferrarese et al., this result for the most luminous AGN indi- 
cates a required total virial mass of 1.5x10^^ Mq for the most luminous obscured galaxies. 
If we assume that 90% of this virial mass is in the form of non-baryonic, dark matter, the 
required stellar mass is 1.5x10^^ Mq, a number only 5 times larger than the median lower 



limit for the observed stellar masses in lBussmann et al.l (120091 ). Starbursts at the maximum 



SFR of 3500 Mq yr ^could accumulate this stellar mass in 5x10'' yr. 

Our various conclusions regarding the luminosities and SFRs of the high redshift AGN 
and starbursts discovered by Spitzer are consistent, therefore, with the interpretation that 
these sources represent the most massive galaxies in the process of formation. Results are 
also consistent with the interpretation that the central, compact massive objects which power 
AGN were in place by the epoch when these galaxies had maximum luminosity. 
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4. Summary and Conclusions 

Continuum luminosities z/L,^(8 fim) arising from dust emission are given for 115 obscured 
AGN and 60 unobscured (type 1) AGN with 0.005 < z < 3.2 using measurements with the 
Spitzer IRS. Obscured sources are defined by having optical depth > 0.7 in the 9.7 /xm 
silicate absorption feature and unobscured (type 1) AGN by having silicate in emission. 

These AGN include the most luminous sources known in the Universe as measured by 
z/Lj^(8 /im). Luminosity evolution for the most luminous obscured AGN is found to scale as 
(1+z)^^ to z = 2.8. For the more limited sample of unobscured, type 1 AGN, the scahng with 
redshift is similar, but the uL^lS /im) is approximately 3 times more luminous at all redshifts 
compared to the obscured AGN. The evolution factor for AGN is the same as previously 
found in a similar way with IRS results for the most luminous starbursts in the Universe. 

Using both obscured and unobscured AGN having total IRAS fluxes, empirical cahbra- 
tions are found between uL^lS /im) and Lir. Combining these calibrations with the redshift 
trends and using a similar result previously determined for starbursts, we conclude that total 
infrared luminosities for the most luminous galaxies scale with redshift as: 

log LiniAGNobscured) = 12.3±0.25 + 2.6(±0.3)log(l+z) for L,^ in L^, 

log L/r(AGN1) = 12.6±0.15 + 2.6(±0.3)log(l+z), and 

log L/ij(starburst) = 11.8±0.3 + 2.5(±0.3)log(l+z). 

Comparisons of obscured and unobscured AGN are consistent with both types having the 
same total luminosities with differences arising only from orientation effects, as expected 
in the unifled model for AGN. The obscured AGN contain very dusty clouds which extinct 
about 50% of the intrinsic luminosity at 8 fj,m. 

The luminosity function of obscured AGN is determined for the highest luminosities 
at z = 2.5, and this luminosity function agrees with that for local IRAS ULIRGs if they 
are scaled by (1+z)^-^ luminosity evolution. At the highest luminosities (Lir > lO^^Lo), 
obscured and unobscured AGN together contribute about 10 times more luminosity density 
to the early Universe than do the most luminous starbursts. 

For comparison with source fluxes found in wide-fleld Spitzer MIPS surveys, extrapola- 
tions of observable fi,{24 /im) to z = 6 are made using the most luminous examples currently 
known. Both obscured and unobscured AGN should be detected to z ~ 6 with f^(24 /im) > 
0.3 mJy, even without continued luminosity evolution for z > 2.5. Type 1 AGN discovered 
using dust luminosity observed at 24 fim show a maximum in luminosity at z ~ 2.5, similar 
to the redshift of maximum luminosity derived from optically-discovered quasars. In contrast 
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to AGN, the most luminous starbursts cannot be detected with fi/(24 fxm) > 0.3 mJy for z 
> 3, even with continued luminosity evolution. 

The maximum luminosities for dusty AGN at z = 2.5 imply a minimum accreting mass 
of 3x10^ Mq. The observed stellar masses of such galaxies could accumulate in < 5x10^ yr 
at the maximum star formation rates which are observed for dusty starburst galaxies. 

We thank L. Sargsyan for help with spectral analysis and SDSS spectra, D. Farrah 
for providing a list of SWIRE Lockman sources, and D. Barry for continuing help with 
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of Technology, under NASA contract 1407. Support for this work by the IRS GTO team 
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Table 1. Continuum Luminosities for Obscured AGN 



Source 


Name 


J2000 coordinates 


z 


f,(8Mm)* log 


; iyU{8 fim )^ 


Ref.^ 


1 


IRAS 00091-0738 


001143.3-072207 


0.118 


52 


44, 


.79 


2 


2 


IRAS 00188-0856 


002126.5-083926 


0.128 


53 


44, 


.87 


2 


3 


IRAS 00397-1312 


004215.5-125602 


0.262 


155 


45, 


.99 


1 


4 


IRAS 01004-2237 


010250.0-222157 


0.118 


92 


45, 


.04 


2 


5 


IRAS 01166-0844 


011907.7-082910 


0.118 


23 


44, 


.44 


2 


6 


IRAS 01298-0744 


013221.4-072908 


0.136 


83 


45, 


.12 


2 


7 


IRAS 06035-7102 


060254.0-710310 


0.079 


132 


44, 


.84 


1 


8 


IRAS 06206-6315 


062101.2-631723 


0.092 


56 


44, 


.60 


1 


9 


IRAS 08572-F3915 


090025.4-F390354 


0.058 


560 


45, 


.20 


2 


10 


UGC 5101 


093551.7+612111.3 


0.039 


278 


44, 


.54 


1 


11 


FSC 09425-f-1751 


094521.4+173753 


0.130 


47.9 


44, 


.85 


3 


12 


NGC3281 


103152.1-345113 


0.0107 


1900 


44, 


.26 


4 


13 


SWIRE 


103303.37+572050.8 


0.50 


2.5 


44, 


.78 


6 


14 


IRAS 10378-^1108 


104029.2+105318 


0.136 


28 


44, 


.65 


2 


15 


SWIRE 


104314.93+585606.3 


2.25 


1.3 


45, 


.78 


12 


16 


SWIRE 


104325.62+581852.3 


0.67 


3.5 


45, 


.19 


6 


17 


SWIRE 


104407.97+584437.0 


0.56 


3.5 


45, 


.03 


12 


18 


SWIRE 


104528.29+591326.7 


2.31 


3.4 


46, 


,21 


12 


19 


SWIRE 


104605.56+583742.4 


2.12 


1.2 


45, 


.70 


13 


20 


SWIRE 


104627.26+592526.9 


0.67 


1.6 


44, 


.85 


6 


21 


SWIRE 


104652.52+571154.9 


1.0 


4.9 


45, 


.69 


6 


22 


SWIRE 


104847.15+572337.6 


1.47 


4.5 


45, 


.98 


8 


23 


SWIRE 


105359.22+592103.0 


0.95 


2.6 


45, 


.37 


6 


24 


SWIRE 


105404.32+563845.6 


1.72 


6.5 


46, 


.27 


8 


25 


IRAS 11095-0238 


111203.4+020422 


0.106 


83 


44, 


.90 


2 


26 


IRAS 11130-2700 


111531.6-271623 


0.136 


25 


44, 


.60 


2 


27 


FSC 11257-1-5113 


112832.7+505721 


0.197 


10.5 


44, 


.56 


3 


28 


IRAS 12018-M941 


120424.5+192510 


0.169 


43 


45, 


.03 


1 


29 


IRAS 12071-0444 


120945.1-050113 


0.128 


66 


44, 


.97 


1 


30 


IRAS 12127-1412 


121518.9-142945 


0.133 


72 


45, 


.04 


2 


31 


NGC4388 


122546.7+123944 


0.0084 


2300 


44, 


.13 


4 


32 


IRAS 12514-^1027 


125400.8+101112 


0.319 


64 


45, 


.78 


1 


33 


Mrk 231 


125614.2+565225.2 


0.042 


1380 


45, 


.31 


1 


34 


FSC 13297-^4907 


133150.5+485150 


0.128 


11.5 


44, 


,21 


3 


35 


Mrk 273 


134442.1+555312.7 


0.038 


287 


44, 


.53 


1 


36 


IRAS 140704-0525 


140931.3+051131 


0.264 


29 


45, 


.26 


1 


37 


NGC5506 


141314.8-031227 


0.0062 


2000 


43, 


.79 


4 
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Table 1 — Continued 



Source 


Name 


J2000 coordinates 


z 


f.(8 Mm)^ 


log vLviS jim )^ 


Ref.'^ 


38 


SST24 


142648.90+332927.2 


1.82 


2.6 


45.91 


13 


39 


SST24 


142731.37+324434.6 


0.90 


1.8 


45.16 


7 


40 


SST24 


142924.83+353320.3 


2.73 


1.5 


45.98 


10 


41 


SST24 


142958.33+322615.4 


2.64 


1.3 


46.02 


10 


42 


SST24 


143001.91+334538.4 


2.46 


5.5 


46.47 


13 


43 


SST24 


143050.84+344848.5 


1.21 


5.0 


45.86 


9 


44 


SST24 


143205.63+325835.2 


0.48 


5.7 


45.10 


5 


45 


SST24 


143234.92+333637.5 


1.12 


5.0 


45.80 


7 


46 


SST24 


143301.49+342341.5 


2.10 


2.0 


45.95 


7 


47 


SST24 


143341.90+330136.9 


0.81 


2.7 


45.25 


9 


48 


SST24 


143508.49+334739.8 


2.00 


3.2 


46.08 


13 


49 


SST24 


143520.75+340418.2 


2.08 


2.2 


45.95 


10 


50 


SST24 


143523.99+330706.8 


2.59 


1.6 


45.97 


10 


51 


SST24 


143539.34+334159.1 


2.52 


3.6 


46.31 


10 


52 


SST24 


143644.22+350627.4 


1.77 


3.0 


45.95 


13 


53 


SST24 


143830.62+344412.0 


0.94 


3.6 


45.50 


9 


54 


IRAS 14378-3651 


144059.0-370432 


0.068 


85 


44.69 


1 


55 


FSC 14448-0141 


144727.5-015330 


0.210 


11.7 


44.66 


3 


56 


FSC 14475-M418 


144954.9+140610 


0.251 


34.8 


45.30 


3 


57 


FSC 14503+6006 


145135.0+595437 


0.577 


22.4 


45.86 


3 


58 


FSC 15065+3852 


150825.4+384122 


0.355 


13.3 


45.19 


3 


59 


IRAS 15225+2350 


152443.9+234011 


0.139 


43 


44.86 


2 


60 


Arp 220 


153457.1+233011 


0.018 


500 


44.13 


1 


61 


FSC 15492+3454 


155108.9+344533 


0.311 


13.6 


45.08 


3 


62 


FSC 15496+0331 


155206.2+032244 


0.193 


28.9 


44.98 


3 


63 


SWIRE 


160532.69+535226.4 


2.75 


4.1 


46.43 


13 


64 


FSC 16073+0209 


160949.7+020130 


0.223 


38.7 


45.24 


3 


65 


IRAS 16090-0139 


161140.5-014705 


0.134 


74 


45.06 


1 


66 


FSC 16156+0146 


161809.4+013922 


0.133 


72.8 


45.05 


3 


67 


FSC 16242+2218 


162626.0+221145 


0.157 


6.4 


44.14 


3 


68 


SWIRE 


164216.93+410127.8 


2.40 


4.2 


46.34 


13 


69 


IRAS 17044+6720 


170428.4+671623 


0.135 


53 


44.92 


2 


70 


SST24 


170916.91+585220.5 


0.24 


1.8 


43.97 


6 


71 


SST24 


170939.20+592728.3 


2.52 


2.3 


46.11 


11 


72 


SST24 


171057.45+600745.2 


2.40 


2.9 


46.18 


11 


73 


SST24 


171115.23+594907.1 


0.59 


9.9 


45.53 


6 


74 


SST24 


171144.26+585225.8 


1.22 


1.4 


45.32 


11 
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Table 1 — Continued 



Source 


Name 


J2000 coordinates 


z 


f.(8 Mm)^ 


log fLy(8 fim )^ 


Ref.^ 


75 


SST24 


171302.37+593611.0 


0.684 


7.4 


45.53 


5 


76 


AMS02 


171315.88+600234.2 


1.8 


0.55 


45.23 


15 


77 


AMS04 


171340.62+594917.1 


1.78 


1.5 


45.66 


15 


78 


AMS06 


171343.91+595714.6 


1.8 


1.9 


45.77 


15 


79 


AMS08 


171429.67+593233.5 


1.98 


1.6 


45.77 


15 


80 


SST24 


171433.17+593911.2 


2.45 


1.95 


46.02 


14 


81 


SST24 


171433.68+592119.3 


1.00 


0.7 


44.84 


14 


82 


AMS09 


171434.87+585646.4 


2.1 


1.0 


45.61 


15 


83 


SST24 


171439.57+585632.1 


1.85 


1.35 


45.64 


14 


84 


SST24 


171510.28+600955.2 


2.60 


2.8 


46.22 


14 


85 


SST24 


171535.78+602825.5 


2.48 


1.5 


45.91 


14 


86 


SST24 


171536.34+593614.8 


2.70 


1.8 


46.06 


14 


87 


SST24 


171538.18+592540.1 


2.65 


3.7 


46.36 


14 


88 


SST24 


171543.54+583531.2 


2.47 


2.6 


46.15 


14 


89 


SST24 


171611.81+591213.3 


2.15 


1.4 


45.78 


14 


90 


SST24 


171750.65+584745.3 


2.43 


6.2 


46.52 


6 


91 


SST24 


171758.44+592816.8 


1.95 


5.1 


46.26 


14 


92 


SST24 


171826.67+584242.1 


1.77 


2.5 


45.87 


14 


93 


SST24 


171844.38+592000.5 


2.08 


4.5 


46.26 


14 


94 


SST24 


171844.77+600115.9 


2.04 


1.5 


45.76 


14 


95 


AMS14 


171845.47+585122.5 


1.79 


1.3 


45.60 


15 


96 


SST24 


171848.80+585115.1 


2.34 


2.3 


46.06 


14 


97 


SST24 


171852.71+591432.0 


0.329 


12.7 


45.11 


5 


98 


AMS15 


171856.93+590325.0 


2.1 


0.6 


45.39 


15 


99 


SST24 


172045.17+585221.4 


3.20 


3.2 


46.43 


14 


100 


AMS18 


172046.32+600229.6 


1.6 


2.0 


45.70 


15 


101 


SST24 


172047.47+590815.1 


2.59 


1.1 


45.81 


14 


102 


SST24 


172048.02+594320.6 


2.29 


2.2 


46.02 


11 


103 


SST24 


172051.48+600149.1 


1.97 


1.3 


45.68 


14 


104 


SST24 


172100.39+585931.0 


1.84 


1.6 


45.71 


14 


105 


SST24 


172119.46+595817.2 


1.64 


3.0 


45.89 


14 


106 


SST24 


172126.42+601646.1 


1.90 


1.5 


45.71 


14 


107 


SST24 


172422.10+593150.8 


2.13 


1.8 


45.88 


14 


108 


SST24 


172428.44+601533.2 


2.43 


1.65 


45.94 


14 


109 


SST24 


172448.65+601439.1 


2.40 


3.0 


46.19 


11 


110 


FSC 17233+3712 


172507.40+370932.1 


0.702 


8.1 


45.60 


3 


111 


SST24 


172542.34+595317.5 


0.43 


2.6 


44.66 


6 
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Table 1 — Continued 



Source 


Name 


J2000 eoordinates 


z 




log i^L„{8 fim )^ 


Ref.'= 


112 


IRAS 20100-4156 


201329.5-414734 


0.130 


77 


45.05 


1 


113 


IRAS 20551-4250 


205826.8-423900 


0.043 


269 


44.62 


1 


114 


NGC7172 


220201.9-315211 


0.0087 


2500 


44.19 


4 


115 


IRAS 23498+2423 


235226.0-1-244016 


0.212 


33 


45.13 


1 



^Observed flux density in mJy at 8 /im continuum peak measured from published spectra or 
from our own extractions, as described in section 2.3. Typical uncertainties are ± 5% for sources 
with fi/(8 ^m) > 5 mJy and ± 10% for sources with fi/(8 fim) ~ 1 mJy. 

''Rest frame luminosity J^Ly(8 fim) in ergs determined using peak f,y(8 /im) and luminosity 
distances from E.L. Wright, http://www.astro.ucla.edu/~wright/CosmoCalc.html, for Hq = 71 
kms-iMpc"\ f7ji/=0.27 and nA=0.73. (Log [uL^{8 /zm)(LQ)] = log [uL^{8 /zm)(ergs s"^)] - 
33.59.) 

Sargsvan et al. I (l2008h : 4 



1 = Farrah et al. (2007); 2 



Shi et al 



Imanishi et al.l (|2007D : 3 
Weedman and Houckr(l2009h : 6 = new spectra in Figure 2 [sources 73,111 from archival 



(|2006l) : 

program 30447 (G. Fazio, P.I.), source 90 from archival program 20629 (L. Yan, P.I.), source 16 



tions, program 50031]; 7 =lBrand et al 



10 



Polletta et al 
(l2008h . 



Si) ; 



14 



Houck et all (120051) : 11 =IWeedman et al.l (l2006al) ; 12 



Yan et al.l (|2007f l: ISaiina et al.l (|2007[) ; 15 



2009f); 9 = Brand et al. 


( 


2008b) 


Weedman et al. 


(2006c 


) ; 13 = 


: 15 = Martinez-Sansigre et al. 
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Table 2. Continuum Luminosities for Silicate Emission AGN 



Durce 


Name 


J2000 coordinates 


z 


f,(8Mm)^ log 


z/L^(8 iJ.m )^ 


log Lir" 


Ref. 


1 


I Zw 1 


005334.94+124136.2 


0.061 


270 


44.93 




16 


2 


PG0157+001 


015950.21+002340.6 


0.1630 


36 


44.93 


46.38 


17 


3 


SWIRE 


021640.72-044405.1 


0.870'' 


9.5 


45.86 




20 


4 


SWIRE 


021729.06-041937.8 


1.146^= 


4.6 


45.78 




20 


5 


SWIRE 


021743.04-043625.1 


0.784f 


3.0 


45.26 




20 


6 


SWIRE 


021808.22-045845.3 


0.712 


5.2 


45.42 




20 


7 


SWIRE 


021830.57-045622.9 


1.40r 


6.2 


46.08 




20 


8 


SWIRE 


021938.70-032508.2 


0.435^= 


2.1 


44.58 




20 


9 


SWIRE 


022012.21-034111.8'^ 


0.166 


2.7 


43.82 




20 


10 


SWIRE 


022431.58-052818.8 


2.068*= 


9.7 


46.59 




20 


11 


IRAS07598+6508 


080433.1+645948.6 


0.148 


200 


45.58 


46.14 


1 


12 


PG 0804+761 


081058.60+760242.0 


0.100 


85 


44.86 




16 


13 


PG0838+770 


084445.26+765309.5 


0.1310 


13.6 


44.31 


45.16 


17 


14 


PGlOOl+054 


100420.13+051300.4 


0.1605 


18.8 


44.64 


45.23 


17 


15 


SWIRE 


103724.79+580513.5 


1.51598 


6.1 


46.14 




20 


16 


SWIRE 


103803.38+572701.4 


1.285S 


11.2 


46.26 




19 


17 


SWIRE 


104255.68+575550.2 


1.4684S 


5.0 


46.02 




20 


18 


SWIRE 


104705.12+590728.5 


0.3913S 


3.2 


44.66 




20 


19 


SWIRE 


105106.20+591625.1 


0.7676S 


4.0 


45.37 




20 


20 


SWIRE 


105158.57+590651.5 


1.81318 


5.7 


46.25 




20 


21 


SWIRE 


105404.14+574019.8 


1.1021S 


5.2 


45.80 




20 


22 


SWIRE 


105959.96+574847.8 


0.4529S 


3.7 


44.86 




20 


23 


SWIRE 


110223.65+574435.9 


0.22628 


3.3 


44.18 




20 


24 


PG1126-041 


112916.66-042407.6 


0.0600 


41.7 


44.11 


44.96 


17 


25 


PG 1211+143 


121417.70+140312.6 


0.081 


106 


44.77 




16 


26 


3C 273 


122906.70+020308.6 


0.158 


220 


45.69 


46.41 


16 


27 


PG1229+204 


123203.60+200929.2 


0.0630 


31.6 


44.03 


44.96 


17 


28 


PG1244+026 


124635.25+022208.8 


0.0482 


17.7 


43.55 


44.50 


17 


29 


PG1302-102 


130533.01-103319.4 


0.2784 


23.9 


45.24 


46.17 


17 


30 


PG1351+640 


135315.83+634545.6 


0.088 


66 


44.64 




16 


31 


PG1411+442 


141348.33+440013.9 


0.0896 


68.8 


44.68 


44.96 


17 


32 


SST24 


142614.87+350616.5 


0.217 


3.0 


44.10 




5 


33 


PG1426+015 


142906.59+011706.5 


0.0865 


61.7 


44.60 


45.16 


17 


34 


SST24 


143132.17+341417.9 


1.037 


7.0 


45.88 




5 


35 


SST24 


143156.40+325138.1 


0.412 


3.1 


44.70 




5 


36 


SST24 


143157.96+341650.1 


0.715 


7.8 


45.60 




5 


37 


SST24 


143310.33+334604.5 


2.4 


12.7 


46.82 




5 
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Table 2 — Continued 



Source Name J2000 coordinates z f,y(8 /im)^ log i^Liy(8 ^im)^ log L^^^^ Ref."^ 



38 


SST24 


143409. 54H 


h334649.4 


0.216 


3.1 


44 


11 






5 


39 


PG1440+356 


144207.46H 


h352622.9 


0.0791 


61.7 


44 


52 


45 


18 


17 


40 


SWIRE 


155936. 14H 


h544203.6 


0.3077S 


4.4 


44 


59 






20 


41 


SWIRE 


160128.58H 


h544521.3 


0.7277*5 


6.9 


45 


56 






20 


42 


SWIRE 


160655.38+534016.9 


0.2137s 


3.9 


44 


20 






20 


43 


PG1613+658 


161357.18+654309.6 


0.1290 


62.4 


44 


96 


45 


65 


17 


44 


SWIRE 


163111.31+404805.2 


0.2575s 


6.8 


44 


61 






20 


45 


PG1700+518 


170124.80H 


h514920.0 


0.2920 


57.6 


45 


66 


46 


37 


17 


46 


SST24 


171124.25H 


h593121.7 


1.49048 


4.2 


45 


96 






20 


47 


SST24 


171233.48H 


h583610.5 


1.663"= 


3.4 


45 


96 






20 


48 


SST24 


171313. 96H 


h603146.6 


0.105 


5.1 


43 


68 






5 


49 


SST24 


171352.4H 


h584201.2 


0.52P 


9.6 


45 


40 






5 


50 


SST24 


171430.70H 


h584225.0 


0.562^ 


3.6 


45 


04 






20 


51 


SST24 


171708. 67H 


h591341.1 


0.645° 


2.8 


45 


06 






19 


52 


SST24 


171747.58H 


h593258.0 


0.2477s 


2.9 


44 


21 






20 


53 


SST24 


171839. 73H 


h593359.6 


0.382"= 


4.1 


44 


75 






5 


54 


SST24 


171902. 29H 


h593715.9 


0.178"= 


9.1 


44 


40 






5 


55 


SST24 


172238. 73H 


h585107.0 


1.617S 


5.4 


46 


14 






20 


56 


SST24 


172328.40H 


h592947.3 


1.34° 


4.8 


45 


93 






19 


57 


SST24 


172619.84+601600.2 


0.9258 


3.7 


45 


50 






20 


58 


SST24 


172704.67+593736.6 


1.13 


14.7 


46 


27 






5 


59 


PG2112+059 


211452.57+060742.5 


0.466 


40 


45 


92 






18 


60 


PG2349-014 


235156.12-010913.3 


0.1740 


23.7 


44 


81 


45 


61 


17 



'^Observed flux density in mJy in continuum at rest frame 8 /xm measured from publish ed spectra 
in references given. Typical uncertainties are ± 5%. For sources in Schweitzer et al. ( 20081 ) (reference 
17), U{8 A'm) is determined by scaling from pubhshed fi,(6 fj,m) using the average spectra in Figure 2 of 
Schweitzer et al. for sources without PAH (characteristic of the highest luminosity sources), such that 
f^(8 Aini)/f^(6 ^m) 1.12. 

'^Rest frame luminosity i'L^{8 /im) in ergs s"""^ determined using luminosity distances from E.L. 
Wright, http://www. astro. ucla.edu/'^wright/CosmoCalc. html, for Hq = 71 kms~^ Mpc~^, 17m=0.27 
and f2A=0.73. (Log [i^K{8 ^lm){LQ)] = log [vL^{8 Aim)(ergs s'^)] - 33.59.) 



°For sources in Schweitzer et al. ( 2008h (reference 17), L^^ is determined as in Sanders and Mirabel 
( 19961 ) for sources having measured IRAS flux densities at all IRAS wavelengths; L^^ — 2.14x10^® 
DL2[i3.5(f^(12 /im)) + 5.2(f^(25 ^im)) + 2.58(fi.(60 fim)) +(f^(100 fim))] for flux densities in mJy 
and luminosities in ergs s" -1 (Lq= 3.83x10^3 ergs s"^). 
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'^References fo r IRS spectrum: 1 
Hao et alJ |2005l ): 17 



Farrah et all (120071): 5 = IWeedman and Houckl (|2009[ ): 16 = 
Schweitzer et alJ |2008h : 18 = iMarkwic k-Kemper et alJ (|2007h : 19 = new ex- 



tractions of sources in the FLS and SWIRE Lockman fields liaving f,/(24 fj,m) > 5 mJy and i? > 20 
mag. which have optical redshifts and classifications cited as type 1 AGN in NED; 20 = new extractions 
of sources available within Spitzer Legacy Program 40539 (G. Helou, P.I.)i ^ limited sample with 
f^(24 /zm) > 5 mJy, which have optical redshifts and classifications as type 1 AGN cited in NED or as 
determined from our examination of SDSS spectra, not including 2 sources also in reference 19. 



''Optical redshift and classification in lLacv et al.l (|2007( ). 
^Optical redshift and classification in lSimpson et al.l (|2006( ) 



^Optical redshift and classification from SDSS. 
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rest wavelength {/j^rn) 



Fig. 1. — Comparisons of the brightest sources yet discovered by Spitzer and measured with 
the IRS at z ~ 2.5. Spectra are averages for the most luminous type 1 AGN (top), the most 
luminous obscured AGN (middle), and most luminous starburst (bottom). All spectra are 
shown at the actual flux densities measured with the IRS at rest frame wavelengths for z = 
2.5, although spectra are averages of brighter sources to improve displayed S/N, as described 
in section 2.1.1. 
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rest wavelength {/j.m) 



Fig. 2. — Rest-frame spectra of all new absorbed sources in Table 1 which are taken from 
the SWIRE Lockman and First Look Survey fields because they have i? > 20 and f,y(24 /xm) 
> 5 mJy. Of the 47 available IRS spectra of sources with these criteria, 11 show silicate 
absorption with optical depth of the 9.7 /im feature r > 0.7 and weak PAH emission, 
satisfying our definition of an absorbed AGN. Nine of these spectra are shown above with 
source numbers from Table 1; sources 73, 111 are from archival program 30447 (G. Fazio, 
P.I.), source 90 from archival program 20629 (L. Yan, P.I.), source 16 from archival program 
40539 (G. Helou, P.I.), and sources 13, 20, 21, 23, 70 from ou r own new observations, pro gram 
50031; The t wo previously publishe d spectra are source 75 (IWeedman and Houckll2009l ) and 
source 17 in (IWeedman et al.ll2006d ). All spectra are normalized to fi/(8 ^m) = 1.0 mJy, but 
zero points are displaced for illustration. The zero flux level for each spectrum is shown by 
the short horizontal line; scaling of actual flux densities of sources can be determined using 
observed values of f,^(8 /im) given in Table 1. 
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Fig. 3. — AGN continuum luminosity z/L,^(8 ^m) in ergs s~^ compared to redsliift for all 
sources in Table 1 (obscured AGN). All luminosities are measured from sources observed 
spectroscopically with the IRS. Triangles are IRAS ULIRGs with silicate absorption r > 0.7; 
crosses are other sihcate absorbed AGN with r > 0.7. Sohd hne is luminosity evolution to z = 
2.5 by factor determined from most luminous absorbed AGN, excluding the 3 most 

luminous ULIRGs shown as diamonds with squares (the 3 sources with spectra in Figure 5). 
Luminosity uncertainties for ULIRGS arise primarily from flux calibration uncertainties and 
are typically ± 5%; uncertainties for other sources, which have been discovered in Spitzer 
MIPS surveys at 24 /im , arise primarily from low signal to noise in observed spectra and 
are typically ± 10%, which is smaller than plotted symbols. (Log [uL^{8 //m)(L0)] = log 
[uU{8 //m)(ergs s"^)] - 33.59.) 
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Fig. 4. — AGN continuum luminosity z/Li^(8 /im) in ergs s~^ compared to redshift for all 
sources in Table 2 (unobscured, silicate emission AGN). All luminosities are measured from 
sources observed spectroscopically with the IRS. Solid line is luminosity evolution to z = 
2.5 by factor (1+z)^'^ determined in Figure 3 for the most luminous o bscured AGN . Long- 



dashed line is form of evolution for unobscured AGN determined by iBrown et al.l (l2006l ) 



using a complete sample of Bootes sources with f;y(24 fim) > 1 mJy and optical redshifts 
and classifications as type 1; line is normalized to the brightest of these AGN observed with 
the IRS. Luminosity uncertainties arise primarily from low signal to noise in observed spectra 
and are typically ± 10%, which is smaller than the plotted symbols. (Log [z/L,^(8 /im)(L0)] 
= log [iyL„{8 /im)(ergs s^-*^)] - 33.59.) 
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rest wavelength (^m) 



Fig. 5. — Rest-frame spectra of the 3 most luminous ULIRGs in Table 1 and Figure 3 as 
determined by uL^i^S jim). Sources are labeled by the numbers in Table 1. All spectra are 
normalized to f,^(8 /um) = 1.0 Jy, but zero points are displaced for illustration. The zero flux 
level for each spectrum is shown by the short horizontal line; scaling of actual flux densities 
of sources can be determined using observed values of f^{8 /im) given in Table 1. 
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Fig. 6. — Ratio of total infrared luminosity Lj/j to z/Lj,(8 fim) as function of uL^{8 /im) for 
absorbed AGN in Table 1 which have total infrared luminosities measured from IRAS in 



Farrah et al.l (120071 ) or llmanishi et al.l (120071 ) . The two most luminous sources are sources 



3 and 32 in Figure 5. Median ratio of 0.95±0.25 indicates that the typical correction for 
absorbed AGN is log Ljji = log lyL^i^S /im) + 0.95±0.25 in ergs s~^, or log L//j = log uL^{8 
/xm) - 32.63±0.25, for L/^ in Lq and z^Lj^(8 /im) in ergs s~^ as used in Figure 3. Error bar 
is uncertainty in ratio produced by observational uncertainty in uL^lS //m). 
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Fig. 7. — Ratio of total infrared luminosity L//j to //Lj,(8 yum) as function of uL,^{8 fim) for 
silicate emission, type 1 AGN in Table 2 which have total infrared luminosities measured at 
all IRAS wavelengths; source luminosities are in Table 2. The median ratio Lir/uL^^S fim) 
— 0.74±0.15 indicates that the typical correction for unabsorbed type 1 AGN is log L/r = 
log i/L^(8 ^m) + 0.74±0.15 in ergs s~\ or log Lir = log iyL^{8 /xm) - 32.84±0.15 for Ljr in 
Lq and z^Li,(8 fj,m) in ergs s~^ as used in Figure 4. Asterisks show values for sources without 
measured PAH features and diamonds for sources with weak PAH features. Error bar is 
uncertainty in ratio produced by observational uncertainty in i'L^{8 /im). 
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Fig. 8. — Redshifts and luminosities for highest redshift obscured AGN. Luminosities are 
LjR in Lq using 1^1,^(8 /xm) from Table 1 and transformation log Lm = log z/Lj^(8 fim) 
- 32.63±0.25, for Ljr in L0 and i^L,^(8 //m) in ergs s~^. Space densities and luminosity 
functions are determined for intervals 1.7 < z < 2.2 and 2.3 < z < 2.7 using distributions 
of points shown. Error bar is uncertainty in ratio produced by observational uncertainty in 
uL^lS /im). 
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Fig. 9. — Luminosity function at z = 2.5 determined from obscured AGN in Figure 8. Space 
densities are in number Mpc~^ per interval 0.4 in log L//?, for L//j in L©. Error bars show 
the differences in results for the two independent redshift windows 1.7 < z < 2.2 and 2.3 < 
z < 2.7. Solid line is th e bright end of the local luminosity function for all IRAS galaxies 
from lSoifer et al.l (119871 ). scaled to z = 2.5 using an evolution factor (l+z)^'^. 
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Fig. 10. — Predicted fluxes measured with MIPS at 24 nm as a function of redshift for the 
most himinous sources known at z ~ 2.5 when extrapolated to higher redshifts. Asterisks 
are predictions for the most luminous, unobscured, type 1 AGN; crosses are predictions for 
the most luminous obscured AGN; and diamonds are predictions for the most luminous 
starbursts, all scaled to the brightest examples shown in Figure 1. For each set, the upper 
curve shows expected fi/(24 /im) if source luminosities continue to scale as for 2.5 

< z < 6, and lower curve shows expected fi/(24 fj,m) if there is no luminosity evolution for z 
> 2.5. Horizontal line is MIPS limit of fi/(24 fj,m) — 0.3 mJy for wide area surveys such as 
in Bootes, the FLS, and SWIRE. Squares are type 1 AGN discovered at 24 fj,m in Spitzer 
surveys having optical classifications and redshifts, as cited in text. 



